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Abstract 

We report on density matrix renormalization-group and analytical work on 
5 = 1 antiferromagnetic Heisenberg spin chains. We study the finite size 
behavior within the framework of the non-linear sigma model. We study 
the effect of magnon-magnon interactions on the finite size spectrum and 
on the magnetization curve close to the critical magnetic field, determine 
the magnon scattering length and compare it to the prediction from the 
non-linear a model. 

PACS Numbers: 75.10.-b, 75.10.Jm, 75.40.Mg 
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I. INTRODUCTION AND CONCLUSIONS 



Since Haldane pointed out the difference between integer spin and half integer spin 



of the gap has been accurately computedH'H. With open boundary conditions (OBC), the 
Hamiltonian: 



gives low energy effective S=l/2 spins localized near the chain ends whose interaction 
decreases exponentially with system size, resulting in an exponentially low lying excited 



Many low energy features of the model can be understood using an approximate 
mapping onto the non-linear a model (NLcrM) . In particular, the low lying excited states 
are constructed from a triplet of interacting bosons (magnons) . There are various types of 
predictions from the NLcrM for the behavior of S=l chains that have been tested to varying 
extents numerically or experimentally. Some of the predictions are not really specific to the 
NLcrM but would follow from any relativistic quantum field theory with a triplet of massive 
particles. These include the relativistic dispersion relation^, exponential temperature (T) 
dependence of the specific heat and Bose condensation transition at a critical magnetic 
field0lll. A qualitative prediction specific to the NLcrM is the absence of boundstates. 
Quantitative predictions, following from the details of the magnon interactions, have 
also been made concerning the twolli and three0 magnon contributions to the neutron 
scattering cross section. However, perhaps mainly because these contributions are so 
small, they have, so far, not been tested numerically or experimentally. 

In this paper we present numerical results on low lying states of the spin chain with 
both open and periodic boundary conditions using the density matrix group (DMRG) 




L-l 



(1.1) 



state.fr^. 
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technique.Ej We keep m=800 states in our DMRG calculations. The biggest truncation 
error is 10~^ and the largest system size is L=200. These results are used to study the 
interactions of the edge spins with each other and with bulk magnons, as well as the 
interactions of 2 magnons with each other. In particular, we determine numerically one 
parameter which characterizes the low energy magnon interactions, the scattering length 
for two parallel spin magnons. We deduce how this single parameter determines finite 
size corrections to the S=2 two magnon state with both periodic and open boundary 
conditions. We also show that the same parameter determines the leading correction to 
the square root singularity of the magnetization at the (lower) critical field and generalize 
this results to finite T and length L. Numerical results on all these quantities give a 
consistent value for the magnon scattering length of the S=l chain, of about -2 lattice 
constants or about — ^/3 where ^ is the correlation length. On the other hand, the exact 
S-matrix for the NLcrM gives a scattering length of —2^ /it, roughly twice as large. Thus, 
the NLcxM does not fair well in its first quantitative comparison to the S=l chain. This 
is perhaps not terribly surprising given that the mapping only becomes exact, at low 
energies, at large S. Better agreement could be expected for larger S or for S=l with 
a (non-frustrating) ferromagnetic next nearest neighbor coupling added to decrease the 
bare coupling constant in the NLcrM. This rather poor agreement suggests that other 
predictions of the NLaM, such as the twolil and three magnonlli contributions to the 
neutron scattering cross section may not be very accurate either. 

The L-dependence of the single magnon energy for open chains has been found to 
behave as 1/L for medium chain length and as A + 0(1/L^) for long chain lengthEl. In 
this paper we will show that those two behaviors for different range of lengths have a 



With periodic boundary conditions (PBC) finite size corrections to the groundstate 
and single magnon energies are exponentially small but, as we show, the lowest 2- magnon 



unified expression W + v'^ sin^ j^. 
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state (with S=2) has energy: 



E2P-Eop = 2^A^ + v^sm^j-^. (1.2) 
We will use the physical quantities obtained for S=l spin chains in previous studies A = 



0.4105 and v = 2.49. Note that the Lorentz invariant dispersion relation E = y + {vky 
has been modified by the replacement: 

k —>■ sin k. (1.3) 

(Throughout this paper we use wave- vectors shifted by tt so that the minimum energy 
single magnon excitation occurs at k = 0.) This modification only affects the expansion 
in 1/L at 0(1/L^). To order 1/L^, this is just the energy of 2 free magnons with wave- 
vectors ±7r/L. The fact that these wave- vectors are non-zero and different reflects the 
hard-core boson approximation in which the wave-function is approximated by a free 
fermion (Bloch) wave-function, multiplied by a sign function to correct the statistics. 
The antisymmetric nature of the Bloch wave-function requires the use of two different 
wave-vectors. Periodic boundary conditions then requires them to be odd multiples of 
±7c/L. a is a new parameter which we introduce, the magnon-magnon scattering length. 
This leads to an 0(1/L^) correction: 

._.„,.,^,i(!^)%(-)(^)%o,i/n ,1.4) 

The replacement of k by sin k is not derived here but is just an empirical improvement 
which reflects the presence of a lattice. In particular, if we ignore A this gives the standard 
spin-wave theory formula. It gives a fairly good fit to the magnon dispersion relation over 
a large range of /£! although it fails badly near k = it (that is k = before the wave- vector 
shift). 

For the case of open boundary conditions (OBC), we derive the following formulas for 
the L-dependence of the energies of the lowest energy states of spin 8=0,1,2,3. 



E, = eo{L - 1) + A, - ^ """"^^^^ 



E, = eo(L - 1) + A, + exp(-^^), 



E3~ Ei = Ja^ + vHm\ ^ +jA^ + vHm^ 



L — 2ab — a V ^ ~ — a 

We use the previously obtained parameters: Cq = —1.401484 and A;, = —0.193166. The 
exponentially small terms in i?o,i arise from the interaction between the effective S=l/2 
edge excitations. E2 — Ei is the single magnon energy. Note that it is approximately 



A^ + {vky with k = 7r/{L — 2ah) up to corrections of 0{1/L^). This is simply the 
energy of a free massive particle in a box with an interaction at the boundaries producing 
a scattering length a^. The formula for E^ — Ei gives the lowest energy two magnon 
state. Note that the second wave-vector occurring here is twice as large. This reflects 
the hard-core boson approximation. The appropriate value of this second wave-vector is 
again determined from consideration of the boundary interactions and contains the same 
boundary scattering length a^. The same inter-magnon scattering length, a, appears as 
in the periodic case. These formulas are only expected to be completely correct up to 
0(1/ L^). The same empirical replacement of k by sin A; has been made. 

The magnetization per unit length, at T = 0, L = 00, close to the critical field is given 

by: 



M/L = — 

TTV 



;i.6) 



(We adopt units where gfiB = 1-) The square root singularity was first proposed by 
Tsvelik,0 using an approximate fermionic representation for the S=l chain. It was later 
arguedlii to be the universal behavior of bosons with repulsive interactions, and there- 
fore to be an exact result for integer spin chains. The linear correction to this formula 
was derived by Okunishi et al.lii recently by assuming a 5-function interaction between 



magnons, c5{xi — Xj), with a replaced by 

a -v^Ac. (1.7) 

We argue here that this formula obtains for any short range interaction (which does not 
produce boundstates), the details of the interaction determining the scattering length. 

The scattering length is defined by the behavior of the (symmetric) scattering phase 
shift in the limit of zero wave-vector. For a symmetric wave-function describing the 
relative motion of 2 particles, the long distance behavior is written in terms of the phase 
shift, 6{k) as: 

ip{x) ^ sm[k\x\ + S{k)]. (1.8) 

As /c — 0, for general short range potentials with no boundstates in the limit A; ^ 0: 

6{k) -ak. (1.9) 

We note that the scattering length can be positive or negative. An infinite hard core 
potential gives a positive scattering length equal to the core size. On the other hand a 
repulsive (5-function potential gives a negative scattering length given by Eq. (\l.7\). 

We also derive the generalization of Eq. ( p..6|) at low temperature T and large L in 



order to fit recent Monte Carlo results.^ We obtain a consistent value of a of about 
—0.34^ from all three fits, as mentioned above. However, using a product wave-function 
renormalization group method to calculate the magnetization, Okunishi et al. obtained 
a considerably larger value, a ~ —0.54^, which is closer the prediction of the NLcxM. 
From the viewpoint of the renormalization group treatment of the one dimensional Bose 
condensation transition,^ we may regard a as the leading irrelevant coupling constant. 
Although its value is non- universal, many different quantities can be expressed in terms of 
it. Thus it plays a similar role to one over the Kondo temperature in the Kondo problem. 
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The boundary scattering length, Ob, resulting from the interaction of bulk magnons 
with the chain end and the effective S=l/2 degree of freedom residing there, is found to 
have a value of approximately -1 lattice constants. 

In the next section we discuss the simpler case of periodic boundary conditions, de- 
riving the wave-functions and energies for dilute bosons with short range interactions in 
terms of the scattering length and obtaining the scattering length for the NLaM. We 
compare Eq. ( |1.2D to DMRG results to determine a. In Sec. 3 we discuss the case of 



open boundary conditions, deriving Eqs. (|1.5| ) and comparing them to DMRG results, 
obtaining at, and a consistent value of a. In Sec. 4 we discuss the magnetization, deriving 



Eq. ( |1.6|) and its finite T and L generalizations, comparing to Monte Carlo results and 
again obtaining a consistent value of a. 



II. DILUTE BOSONS AND THE NON-LINEAR a MODEL 



Following the derivation of the NLcrM for large S from the Heisenberg mo del0'0, 
first define 



we 



S2i-i = -S(f)i + li, S2i = S(j)i + li, i = ltoL/2, 



(2.1) 



for a chain of L sites, (p ^i-nd 1 represent the low energy Fourier modes of the spin operators 
with wave- vectors near vr and respectively. Starting from Eq. (|1 . 1| ) we obtain a continuum 
NLaM Hamiltonian: 

2" 



n 



dx 



il2 

S ^ A\dx 



+ ... 



(2.2) 



This gives a velocity of 25* and a coupling constant g = 2/S. The non-linear constraints 
and commutation relations give the corresponding Lagrangian: 



(2.3) 



with the constraint (j)^ = 1. (We have set the effective "velocity of hght" v to one.) The 
spectrum of this field theory is known to consist of only a triplet of massive bosons, created 
by the fields, 0, and their multi-particle scattering states. There are no boundstates. 
The boson mass is exponentially small in the coupling constant, g\ it results from non- 
perturbative effects. A simpler Lagrangian which is expected to have qualitatively similar 
physics is the 0^ model in which the constraint on (f) is relaxed: 

C = Uj-d^$-^$'--^$\ (2.4) 

This model also has a triplet of massive bosons and presumably no boundstates for A > 
(repuhlusive interactions). Of course, the details of the boson interactions will be 
somewhat different in the two models. However, the single boson energy, /S? + v'^k?, for 
A the renormalized mass, is independent of the interactions and simply reflects Lorentz 
invariance and the stability of the single boson excitation. 

We now consider the low energy states of a dilute gas of bosons, all with S"^ = +1 
such that the average spacing is much greater than the Compton wavelength of the boson, 
^ = v/A. A great deal of universality occurs in this limit. In particular since the bosons in 
these states have small wave-vectors, k << C,~^, the usual non-relativistic approximation 
to the dispersion relation is adequate: 

e(A:)^A + ^. (2.5) 

This implies that the behavior will not depend on the Lorentz invariance of the underlying 
field theory. This is an important point because the S=l chain is only approximately 
Lorentz invariant. Nevertheless, it contains stable bosonic magnons whose dispersion 
relation is given by Eq. ( |2.5| ) for small k. (Note that we have shifted the wave-vectors 
in the spin chain problem by vr.) This follows from the assumption that the dispersion 
relation is an even function of k analytic near k = and defines the parameters A and 
V. Due to spin conservation, the low energy states under consideration are completely 
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stable against decaying into states involving bosons of other spin polarizations, so they 
may be integrated out. At low energies, the effect of all processes involving creation and 
annihilation of virtual particles, as given, for example by the Feynman diagram expansion 
of the 0^ theory, may be boiled down to a simple non-relativistic effective Hamiltonian, 
written in first quantized notation as: 



Here m = A/v"^ (the Einstein relation) and V is some effective short range interaction. 
We emphasize that the exact form of V is unknown in all cases except for the weak 
coupling limit of the 0^ model where it is simply a repulsive S function. We do know 
that it is such as not to produce any boundstates and we expect its range to be roughly 
C,. We may indirectly infer its properties, in the case of the NLaM, from the low energy 
S-matrix which is known exactly. Nothing is known, a priori about the effective potential 
for the S=l chain, but we may indirectly deduce some of its properties from numerical 
simulations. Importantly, in the dilute limit the spectrum only depends very weakly on 
the detailed form of V. 

In this dilute limit we may construct the many body wave-functions following the 
construction of Lieb and Liniger for the 5-function Bose gas.H This follows because we 
may ignore the possibility of more than 2 particles being within a distance ^ of each other 
simultaneously so we only have to consider 2-particle scattering processes. In the most 
likely case where all particles are far from each other compared to ^, we may write the 
N-body wave-function as a sum of plane waves. Considering the case Xi < X2 ■ ■ ■ < xn, 
we write: 



for some set of wave- vectors, kj. Here P permutes the fcj's and the sum is over all per- 
mutations, ip is determined for other orderings of the Xj from the symmetry of the wave- 




(2.6) 



N 




(2.7) 
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function following from Bose statistics. By considering what happens when 2 particles 
approach each other, we can see that 

A{Q) = _^(p)e^25[{fc»-fc.)/2]_ (2.8) 

Here the two permutations P and Q differ by the interchange of only 1 pair, i and j and 
6{k) is the even channel phase shift defined by the 2-body problem with potential V{x). 
The Hamiltonian for this problem, where we consider only the relative motion of a pair 
of particles, defining x = xi — X2, is: 

(Note that the reduced mass, m/2 occurs here.) The even wave-functions, at distances 
>> C, take the form of Eq. (|1.8| ), defining 6{k). In general, the coefficients A{P) in 
Eq. (p.7|) are just products of the factor in Eq. ( p.8|) over all all elementary permutations 
corresponding to P. By considering again, a region where all the particles are far apart, 
so y ~ 0, we see that: 

E^—ye. (2.10) 

3 

Periodic boundary conditions then determine the fc^'s, which must all be different from 
each other, from the N equations: 

(_l)^-ie-'=.^ = exp [iY.^26[{h-k,)/2]^ , all j. (2.11) 

We emphasize that this result only holds in the dilute limit when all the fc^'s are small. 
Thus we are only interested in the phase shift for small k. It is easy to see that 

5{k) -ak, (2.12) 

for some constant a at small k, k « 1/^. This follows from observing that in this limit, 
in the region ^ << |x| << 1/fc, the wave-function obeys approximately the zero energy, 
zero potential Schroedinger equation: 
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^ ^(x) = 0. (2.13) 



The even solutions have the form: 

i! (x\x\- a. (2.14) 

Matching this with sin[A;|x| + S{k)] at large |x| then determines 6 ~ —ak. By analogy 
with the standard definition in three dimensions, we refer to a as the scattering length. 
We discussed it briefiy in Sec. 1. We generally expect it to be of O(^). 

It is perhaps worth remarking here that, with our definition of phase shift, it takes 
the value 7t/2 for all k for V = since the symmetric wave-functions are then coskx. In 
this non-interacting case Eq. ( p.l4| ) is not valid. The fact that 6{k) — > at /c — > for 
any short range repulsive interaction, no matter how weak, represents a sort of infrared 
singularity of one dimensional single particle quantum mechanics. It is basically this 
singularity which is responsible for the fact that, in the dilute limit, the many body wave- 
function for interacting bosons reduces to that of free fermions times an antisymmetric 
sign function.lll A related observation is that 6{k) = implies the kj must all be different 
since otherwise the wave-function of Eqs. ( pTTl) , ( p.8|) would vanish. 

In general, it is non-trivial to solve Eq. ( p.llD , requiring numerical methods. However, 



in the small k limit, where Eq. ( p.l2| ) applies, it is easy. Since the phase shift for the 
5-function potential and the hard core potential have this form at low k, we simply recover 
the low density, low k limit of those problems.^ In that limit, Eq. ( p.ll| ) can be written: 

kj{L — No) + a^^ks = TTUj, allj, (2-15) 

s 

where the integers, rij are all even, for N odd and all odd for N even. In the case N = 2 
the lowest energy solution of Eq. (|2.15| ) is: 



ki = -k2 = (2.16) 
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and the energy is: 



J5«-f— )'«-f-y% — f-V. (2.17) 

m\L-2aJ m\Lj mLKLj ^ ' 



This is fitted to our DRMG results for the excitation energy of the lowest state with 
= 1 with PBC in Fig. |1]. A is clearly discernible and gives an estimate of the 
inter-magnon scattering length: 

a ^ -0.34^ (2.18) 

We emphasize that a negative scattering length does not imply an attractive interaction. 
In particular, as mentioned in Sec. 1, a repulsive (5-function interaction gives a < 0. 

The exact S-matrix has been conjectured for the 0(n) NLcrM based on factorization 
of m-particle scattering and various consistency conditions and checked against the 1/n 
expansion. ii The three possible states of a single magnon are labeled by a vector index 
i = 1, 2, 3 so that \Ai >oc 0*|O >. 0(3) symmetry then implies that the S-matrix takes 
the form: 

ikSji = < Aj{p[)Ai{p2),out\Ai{pi)Ak{p2),in > 

= S{pi -Pi)S{p2 - P2)[SikSjicri{9) + 5ij5ki(J2{9) + 6ii6jkCT:i{9)] + {i ^ k,pi ^ P2). (2.19) 
Here the momenta of the particles are labeled by the rapidities 6a'- 

vpa = Asmhea, (2.20) 
and in Eq. ( p.l9| ) the (Xi's are expressed as functions of 



These functions are given by: 

(T2{9) 



9 = 6,- 02. 




2ni9 




{e + m){e- 


2m) 


9(6 - ni 


) 


{e + m){e- 


27ri) 




ni) 



[2.21) 
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The S-matrix for scattering of 2 spin-up particles is obtained by using: 

\A+{p) >= > +t\A2{p) >)/V2, (2.23) 



giving: 



< A+{p[)A+{p'^),ont\A4p,)A+{p2),m > 

= -Kvi - p'Mp2 - P2)\^^ - 5(P2 - p'Mp, - P2)\i^- (2.24) 



This determines the phase shift: 



,2iS[{pi^P2)/2] 



1 + i9/7i 



(2.25) 



l-i9/7r' 

Now letting pi — p2 ^ 0, we obtain the exact scattering length of the NLcrM, for 2 spin 
up magnons: 

aNLaM = ~2v/An = -2^/7i, (2.26) 

roughly twice as large as the value determined numerically for the S=l Heisenberg model. 

This calculation can also be repeated for 2 magnons in a state of total spin S=l or S=0. 
The S=l states are produced from the groundstate by the low wave- vector components 
of the lattice spin operators. We obtain the scattering lengths for spin-S from the NLaM: 

= -2e/vr 
a' = -e/vr 

a° = e/vr. (2.27) 

For the S=0 case, where the wave-function is again symmetric with respect to the spatial 
co-ordinates, Eqs. (|2.16|) and ( p.l7|) still apply with a replaced by a°. The S=l case is a 



bit different because now the wave-function must be anti-symmetric with respect to the 
spatial co-ordinates since it is anti-symmetric with respect to the spin indices. In this 
case we obtain the wave-vectors for the lowest energy states: 
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ki^O, k2 ^ ±271 /{L-a^), 



(2.28) 



and the energy: 



£;,_£„,«2A+-(-) +-(-) -. (2.29) 



III. OPEN BOUNDARY CONDITIONS 

While a periodic integer spin chain has an excitation gap which depends only weakly 
on L, in the case of OBC there are low lying edge excitations. These correspond to 
effective S=l/2 spins localized at the two ends of the chain which are coupled to the 
bulk degrees of freedom. An effective field theory is given by the NLaM with Neumann 
boundary conditions, 

d(p/dx = 0, {x = 0,L), (3.1) 

coupled to the S=l/2 spins, Si and S2. In general, one expects a coupling to both the 
staggered and uniform magnetization density at the chain ends:i^ 

H = ^bulk + ^sim ■ Si + (-l)V(^) ■ S2] + Xuim ■ Si + l(^) ■ S2]. (3.2) 

Integrating out the bulk field, gives an effective exponentially small interaction, between 
the end spins proportional to A^. 

Je// oc A^e-^/«. (3.3) 

This leads to the splitting between the lowest singlet and triplet states in Eq. (|1.5|) which 
can simply be thought of as formed from the two S=l/2 end excitations. The singlet 
(triplet) state of the two S=l/2's has lowest energy for L even (odd). (We only consider 
L even here.) We check the formulas in Eq. (|1.5| ) for Eq and Ei in Figs. ^ and ^. Because 
Jeff is exponentially small, it becomes smaller than our numerical error for the longest 
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chain lengths that we have studied. Therefore, we do no show data for these chain lengths 
in Figs. ^ and |^ since it would be meaningless. This problem does not occur for the other 
energy differences that we consider which only involve finite size effects scaling as powers 
of 1/L. 

The lowest S = 2 state is obtained by polarizing the two end spins and then adding 
one = 1 magnon. When the end spins are polarized, the term produces an effective 
boundary potential energy acting on the magnons. We may describe this by adding an 
effective boundary potential for the magnons. Thus the non-relativistic N-body Hamilto- 
nian of Eq. (|2.6| ) gets modified to: 

I N j2 N 

Hn = ~^T.^+ H V{x.-x,)+Y.[V,{x,) + V,{L-x,)], (3.4) 

1=1 i <i,j> i=l 

where < Xj < L with the Neumann boundary condition: 

djilj{xi,X2, . . .xn) = 0, {xj=(],L). (3.5) 

We expect Vb{x) to be some short range function (with range of 0(^)). It is not a priori 
clear whether Vb is attractive or repulsive. However, the absence of boundary boundstates 
that is evident from earlier numerical workil indicates that it is repulsive. 

For low energy multi-magnon states only the scattering length a^,, produced by Vf, is 
relevant to the spectrum. This has the effect of modifying the low energy single magnon 
wave-functions, at distances large from the boundaries compared to ^, to the form: 

ip^x) ~ sm[k{x — ab)] oc sm[k{L — ab — x)]. (3.6) 

Consistency of these expressions quantizes the wave-vectors: 

k,=nj/iL-2ab), j = 1,2,3.... (3.7) 

The effective chain length becomes L — 2ab. It is now clear that our results are quite 
robust against variations in how the effective Hamiltonian is written. For instance, it 
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doesn't matter whether we impose the Neumann b.c. at x = or x = 1. Such differences 
can be absorbed into the scattering length. We also note that < Sj > exhibits period 
7T oscillations upon which are superimposed long wavelength variations which can be 
interpreted in terms of the magnon wave function.B The boundary scattering length can 
be determined from the additional energy to add one = 1 magnon to the polarized end 
spins, giving the third of Eqs. (|1.5|) . This gives a term: 



This is equivalent to imposing vanishing boundary conditions on a chain of length L + 2. 
We remark that there is no simple derivation that we know of for this result. Presumably 
the value of depends on the (integer) magnitude of the spin and other details of the 
microscopic Hamiltonian. We note that a somewhat better fit is obtained in Fig. § by 
the replacement k sin k. 

We now consider the spin-polarized two-magnon wave-function in the presence of the 
polarized end spins, with = 3. The effect of the boundary potential is simply to 
fix the effective size of the chain at L — = L' with an effective vanishing boundary 
condition. The 2-magnon wave-function, with wave vectors ki and ^2 can then be written 
in terms of the magnon-magnon phase shift, 6{k). This wave-function is made from linear 
combinations of the periodic 2-magnon wave functions discussed in Sec. 2 with wave- 
vectors {ki,k2) {ki,—k2), {—ki,k2) and {—ki,—k2). Two different magnon-magnon phase 
shifts occur: 




(3.8) 



Fitting to this expression, as shown in Fig. ^, we obtain: 



-1. 



(3.9) 



6± = 6[{k,±k2)/2]. 



(3.10) 



for Xi < X2i this wave-function, constructed to obey ■?/'(0, X2) = 0, is given by: 



16 



1p{Xi X2) = e^^^^^^'^^'^^'^^ — g*[('^2a:i+A:ia::2)+2(5_] _ ^i{-k\x-i_+k2X2) _j_ ^i[{k2X\-k\X2)-25+\ 
j^\^_^i{klx\-k2X2) _|_ gi[(-fc2Xi+A;iX2)+2<5+] _|_ ^-i{kixi+k2X2) _ ^-i[{k2Xi+kiX2)-2S-]-^^2i{5- -5+) -^-^^ 

The wave-function for X2 < xi is obtained by interchanging xi and X2 in order to enforce 
the symmetry required by Bose statistics. Imposing iIj{xi,L') = then requires: 

^ik2L' ^ ^-ik2L'+2i{5--5+) 

Now using the small k approximation to 6{k) gives the conditions: 

2k2{L' - a)k2 = 27^2 

2ki{L' - a)ki = 2nni. (3.13) 
Thus the allowed wave- vectors are: 

h = -r-TT^ — , (3.14) 

L — /aft — a 

The lowest energy state has rii = 1, n2 = 2. Fig. ^ shows that the second magnon has 
momentum k2 = ^ with n2 = 2. [The must be different in order for the wave-function 
of Eq. ( |3.11| ) not to vanish.] Setting ai, = —1, its energy is given by the last of Eq. ( |1.5| ). 
As can be seen from Fig. good agreement with this formula is obtained with a value 
of the magnon-magnon scattering length: 

a ^ -0.32^. (3.15) 

This appears consistent, within the numerical error, with the result a ~ —0.34,^ obtained 
with PBC. 

We also comment briefly on finite size energies of two magnon states with other spin 
quantum numbers. This would certainly be simplest to study using PBC so that there 
are no end spins to worry about. An alternative approach, following White and Huse,i 
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would be to consider OBC but add extra S=l/2 variables at the edges of the system 
so as to cancel the effect of the edge excitations. The boundary scattering length will 
depend on the details of the boundary couplings. The contribution of the magnon-magnon 
interaction energy to the 2 magnon states can again be expressed in terms of the scattering 
length for the appropriate spin channel, . It can be seen that the wave- vectors are still 
given by Eq. ( p.l4|) with the appropriate value of a for both symmetric (S=0,2) and 



anti-symmetric (S=l) wave-functions. Thus the energy, to 0(1/L^) in all cases is given 
by: 







5 






2m 




m 


{£ 


/ L 



Eso-E,„^2A + -{^-) +-[^)^ (3.16) 

where L' = L — 2af, and aj, is the appropriate boundary scattering length. 

DMRG results on 2 magnon energies of various S were reported by White and Huse. 
They found repulsive interaction energy for S=0,2 but attractive for S=l. This is com- 
pletely consistent with our approach. They did not use unequal wave-vectors for the 
2 magnons in their definition of the interaction energy, in the even S case. Thus their 
positive interaction energy is presumably of 0(1/L^) and just reflects the effective Fermi 
statistics of the dilute Bose system, i.e. the requirement of unequal wave-vectors. The 
attractive interaction energy in the S=l channel corresponds to a negative a^. They only 
present data for L=60, but if we assume that the term is dominating at that value 
of L we extract a value, ~ —0.43^, about 30% larger in magnitude than the NLcrM 
prediction. 



IV. MAGNETIZATION 

We now consider wave-functions, with periodic boundary conditions, containing an 
arbitrary number, A^, of magnons, assuming that the density and all wave-vectors are 
small compared to 1/,^. To solve Eq. ( |2.15| ) for the k/s for a general low energy A^-particle 
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state we now use the fact that the density is low, Na/L « 1 and aj^s^s « kjL. Hence, 
in lowest order approximation: 

kj kjo = nrij/L. {rij even for N odd, rij odd for N even) (4.1) 

The dimensionless expansion parameter is na where n is the density, N/ L. Hence we may 
expand the kj in powers of a. The leading correction is kj = kjo + Skj, where 

Sk^ = ^kjo - I ^ kso. (4.2) 

The energy of this state is approximately: 

^ - ^ Uk% + 2%o<^%] - ^ E ^^0 + ^ E (^^0 - %o)^. (4.3) 

3 j <i',3> 

We now consider the groundstate energy in the limit L — > oo with a non-zero but small 
density, n. The values of kjQ occupy the "Fermi sea", \k\ < kp where kp is determined 
from the density: 

f^F dk kp 

n= — = — . (4.4) 
J-kp 2n n ^ ^ 



The groundstate energy, to 0{a) is then: 

^ , /A TT^ f'^" dk k^ a f^F dk f^p dk', 

Eo(n)/L = (A- H)n+ + — / — / — ik-k'f. (4.5) 

^" ^ ' J-k,.2n2m 2m J-kp 2t: J-kp 2t:^ ' ^ ' 

Here we have included the Zeeman term, —Hn in the energy; H is the applied magnetic 
field. Performing the integrals gives: 

E,{n)lL = (A-H)n+^ + (4.6) 

6m 3m 

This represents an expansion of the energy in powers of the density. The first term which 
depends on the interactions is the O(n^) term. Minimizing Eq with respect to n gives the 
density or magnetization per unit length: 



M{H)/L^n^ — 
nv 
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(4.7) 



In Fig.|^ we plot M/L vs. H from Eq. ([4.7|) with and without the interaction (a = —0.34^ 
or a = 0). The leading order correction due to the magnon-magnon interaction is obvious 
around magnetization M/L = 0.02. 

This finite density correction can be generalized to finite temperature, T, as was 
observed by OkunishiS in the special case of the (5-function interaction. The correction 
to the free energy of lowest order in density is given by including thermal occupation 
numbers in the second term of Eq. (|4.5| ). 

AF=—J -^-^nF{h)nFik2)ih - k^f . (4.8) 

Note that it is the Fermi distribution function which appears, rather than the Bose 
function. This simply follows from the condition that the kj should all be distinct so 
that there is an effective occupation number for each momentum which can be or 1 
only, npik) is evaluated at finite T and H, then the magnetization is obtained by the 
usual thermodynamic formula, OF /OH = —M. In order to compare with recent Monte 
Carlo data on the magnetization for the S=l chain, it is useful to also generalize our 
formulae to finite length, L with periodic boundary conditions. There is a slight subtlety 
in doing so because the allowed wave- vectors of the magnons alternate between k = 27m /L 
( "even wave- vectors" ) for an odd number of magnons and k = 27i{n + l/2)/L ( "odd wave- 
vectors") for an even number. This follows from the sign change of the wave-function 
each time one magnon passes another one. However, this is easily dealt with exactly by 
inserting appropriate factors of: 

(-1)^ = e^"S^"^ (4.9) 

into the partition function trace. This effectively gives the chemical potential an imaginary 
part, essentially converting fermion occupation numbers into boson ones. For a = 0, the 
partition function is given by: 

= (1/2)[Z°, - + + 1/ZU (4.10) 
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Here Zp^ denotes the partition function for free fermions with even wave-vectors; 
denotes the partition function for bosons with odd wave-vectors, etc. Note that inverse 
boson partition functions occur. The thermal average of the second term in Eq. ( ^.3| ) 
becomes: 

AF = ± ih - hf[Z%^n,{h)n,{h) 

\ki,k2 

nB{ki)nB{k2)\ 



+ Y.iki-k2f[Zl,nF{ki)nF{k2) 

+^nB{k,)nB{k2)]^ (4.11) 

where the first sum is over even wave-vectors and the second over odd wave-vectors. The 
resulting magnetization curve is plotted in FigJ^. Note the smoothing of the singularity 
at the critical field due to a finite T and the oscillations due to a finite L. We use the 
interaction parameter, a ^ —0.34^, obtained from the ground state numerical data so 
there are no free parameters in drawing Fig.^ The agreement with the Monte Carlo 
data,lii at the same length and temperature, L=100, T=l/100 is remarkable. 

This work is partially supported by Chinese Natural Science Foundation, HKUGC 
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FIG. 1. Excitation energy of two magnons state with PBC, E2P — Eqp. 



Fig.2, Lou et. al., PRB 
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FIG. 2. Ratio of coupling energy of two edge 1/2 spins in ground state and edge excited 
state. Here energy ratio vs. — 1) is plotted. Dashed line is — ^. 
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Fig.3, Lou et. al., PRB 
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FIG. 3. Coupling constant for edge spins. Here ln(i?i — Eq)/[L— 1) vs. — 1) is plotted. 
The straight fitting line is ln(£;i - Eq)/{L - 1) = -l/i + c'/{L - 1), where ^ = 6 and c'=2/3 is 



taken from Ref.[ 16]. 
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Fig.4, Lou et. al., PRB 
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FIG. 4. Single magnon energy with OBC. Here {E2 — Eif' vs. sin^ is plotted. The 
boundary scattering length, 05 ?s — 1 is obtained. The straight fitting line shows a better fit at 
short chain lengths results from the replacement of k by sin k. 
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Fig.5, Lou et. al., PRB 
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FIG. 5. Two magnon excitation behavior for 5 = 1 open chain. Here (£^3 — £^2)^ vs. 
sin^ is plotted. The straight line is {Ez - £2)^ = + sin^ without any behavior of 
magnon-magnon interaction included. 
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FIG. 6. Two magnon excitation energy for OBC. Here — E2 — y A'^ + v'^ sin^ vs. 

l/L^ is plotted. The fittin g curve goes to zero linearly as 0. The coefficient of the cubic 

term determines the scattering length, a ^ —0.32^. 
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Fig.7, Lou et. al., PRB 




FIG. 7. Magnetization curve for S = 1 chain near critical field He = A. (The exchange 
constant and g^B are set equal to 1.) For infinite length and at T = 0, we plot it as given in 
Eq.( |4.7| ). Full line has included the leading order contribution of magnon-magnon interactions. 
Dashed line is a reference line for comparison and it is for free hard core boson approximation. 
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Fig.8, Lou et. al., PRB 




FIG. 8. Magnetization curve for S = 1 chain near critical field He = A, with length 
L = 100 and at temperature kT = I/IOO. Full line has included the leading order contribution 
of magnon-magnon interaction. Dashed line is a reference line for comparison and it is for free 
hard core boson approximation. The dots are the Monte Carlo resultsEl. 



28 



